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t-H Abstract 

Consider {^wj^gx' a branch of unstable solitons of NLS whose linearized operators 
have one pair of simple real eigenvalues in addition to the zero eigenvalue. Under radial 
symmetry and standard assumptions, solutions to initial data from a neighbourhood of 
the branch either converge to a soliton, or exit a larger neighbourhood of the branch 
transversally. The qualitative dynamic near a branch of unstable solitons is irrespective 
of whether blowup eventually occurs, which has practical implications for the description 
p I of blowup of NLS with supercritical nonlinearity. 

< 

1 Introduction 

-t— > 

£h 1.1 Nonlinear Schrodinger equation and solitons 

Let us consider the nonlinear Schrodinger equation (NLS) with covariant general nonlinear- 
itjg(u) :=f(\u\ 2 )u, 

XT) id t u + Au-V (x)u + f{\u\ 2 )u = 0, (1.1) 

where u : (t,i)elxl f 4 u(t,x) £ C, N > 1, and ti(0, •) = u G F 1 ^). We assume 
that / is real-valued, the C 1 nonlinearity g(u) is C 2 for u / 0, and i/^-subcritical in the 
f — sense that: 

o 

5 (0) = 0, 5(0) = and \g"(u)\ < \u\ mi ~ 2 + \u\ m2 ~ 2 for u ^ 0, (1.2) 

with 1 < mi < m% < m max ; m max = if N > 3, and m max = +oo if iV = 1, 2. Vq is 
^ either zero, or a smooth localized potential such that a range of assumptions hold; see 



Section 1.2 These equations are Hamiltonian, with the following conserved quantities: 



E(u ) = E(u) := l -j(\Vu\ 2 + V \u\ 2 )-jG{u) (energy), 



M(uq) = M(u) := - I \u\ 2 (mass), 



where G(s) = Jq g(s)ds. For data in the energy space H 1 ^ 1 ^), equation (1.1) is locally 
well posed with a blowup alternative: there exists T max > such that 

u (t)GC([o,r max ),ii' 1 ) J 
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and either T max = oo or limt_>.T max H^WII/f 1 = 00 • The modern proof is due to Kato |20j. 

Another well known property of is the existence of solitary wave (or soliton) 

solutions of the form u(t,x) = (p UJ (x)e tujt . Berestycki and Lions [I] proved that for Vo = 
there exists a solution (f> u > of 

A&,-w&, + /(|0 w | 2 )0 w = 0, (1.3) 

provided lo > and there exists u\ > such that G{u\) > ^u\. Note that this further 
assumption on the nonlinearity g is independent of N. The positive solution of (1.3) is 
unique and radially symmetric up to translations. See Cazenave's book [8] and the references 
therein, particularly McLeod [28J . For Vq 7^ 0, see Rose and Weinstein [33J and articles that 
cite it. 

We are interested in the dynamics and stability of these objects. Because of rotation 
invariance of equation (1.1), the proper notion is orbital stability: we say that <j) UQ is orbitally 



stable if for all e > 0, there exists 5 > such that 



||«o - ^oWm < # => sup inf \\u(t) - (/) Uo e ts ffl < e. 

Note that, when the soliton family has the spatial translation symmetry, the solution may 
move but remain close to the soliton family. In that case, the above concept of orbital 
stability is not adequate, and the nullspace of the linearized operator (to be discussed in 



Section 2.1) will have extra elements (V0 W ). One possible remedy is to enlarge the "orbit" 
to contain all translated solitons. Another way, which we choose in this paper, is to restrict 
ourselves to the cases with no translation. For example, we may assume the perturbation is 
either radial or even. Alternatively, we may also assume the potential Vq is nonradial and 
hence the soliton family has no translation symmetry. 

We now assume, as is usually the case, that we have a whole branch of solitons, that 
is a family {(j> w ) w gi of solitary waves, with <f> w > radial, X = (w, tJ), and a C 2 map 
lo € X 1— > <\) w £ i7 1 (IR iV ). Under these assumptions and some generic spectral conditions, for 
ujq € X, we have: 

(i) WO is orbitally stable if ± J |0 W | 2 | w=wo > 0; 

(ii) <f) Wo is orbitally unstable if ^ J |</> w | 2 \ UJ=UJo < 0. 

The part (i) is due to Cazenave-Lions [9] and Weinstein |40| . by variational and energy 

methods. The part (ii) is due to Shatah and Strauss [35] . See also [T7]. We are interested 

in unstable branches of solitons. 

In the case of a pure power focusing nonlinearity, f(s) = +s( m ~ 1 )/ 2 and Vq = 0, there ex- 

1 

ists a profile (f> u for all lo > 0. Moreover, there is scaling invariance, 0cj(x) — LJ m 1 <^>i(-y/uJx), 

1 N . 

by which we can compute ||^ w ||x,2 = w™- 1 4 ||0i|| L 2. We let m c = + 1 and note that ^ 
is stable for all to > when 1 < m < m c (I? subcritical case), and (f> u is unstable for all 
lo > when m c < m < m max (L 2 supercritical case). 

When the solitary wave <p uo is stable, it is likely that a nearby solution will relax to 
some (j) u+ as time goes to infinity. The frequency cj+ is close to loq but most likely different. 
In addition, the convergence happens only locally since there is radiation going to infinity 
with unvanishing L 2 -mass. Thus one also considers a local concept of stability of the 
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,j£Z is asymptotically stable if, for any uq in a suitable 



branch: we say that the branch (<f> u 
neighbourhood, 

\\u(t)-^ {t) e te ^\\ L 2 — ►(), as t 

v ' loc 

for some continuous functions uj(t) G X and 0(i) G R. Asymptotic stability is the main 
stability concept used in the statement of our main theorem. 



oo, 



This paper is organized as follows. Section 1.2 will itemize our assumptions in terms of 
the nonlinearity and in terms of the evolution operator JL linearized about cp w . We then 
state our main result. Section 11.31 will discuss the context for our theorem and relevant 



literature. The linearized operator JL will be properly introduced in Section 2.1, followed 
by the proof of our main result. 



1.2 Assumptions & Main Result 

1. Solitons exist for u G I, and, moreover, the profiles 



and their derivatives d^o xYu; 



have exponential decay. We assume the map u i— ^ <p u G H 1 is C 2 . 

For all u el, the linearized operator JL has eigenvalue zero with multiplicity two. 
The discrete spectrum contains exactly two simple real eigenvalues, e+ > > e_, 
with corresponding eigenfunctions 1+ and Y- with exponential decay. The purely 
imaginary continuous spectrum is bounded away from zero. There is no embedded 
eigenvalue or resonance in the continuous spectrum or its endpoints. 

As unstable solitons are typically characterized by {du(j>u,<t>J) < 0, we assume that 
u , 4>uj) 7^ with a uniform bound on I. 



3. Nonlinearity g(u) 



f(\u\ )u is sufficiently strong. In addition to (1.2), we require 
2 



mi > 1 + N (1 + Gp) 



mi > 1 + 



1 + 



N V rn + 1 
where mo = min {2, mi} . 



(1.4) 
(1.5) 



Here, we use p to denote 777,2 + 1 = p < p max = m max + 1, and the standard notation 
a r = N (1 — 1) and ^ + h = 1 for any exponent r. See equation (2.24), below, for a 
related fixed constant q G \p,Pm&x) determined by both mi and mi- The choice of q 
will require N > 2. 

Condition (1.4) implies mi > \ + + ^ n2j ^ n —, the Strauss exponent. However, 
we do not approach this value, as condition (1.5) implies either mi > 1 + % (!) or 



mi > 



iv 



+ 



N 



, which are more restrictive. See figures [l and [2J We note 



that conditions ( |1.4[ ) and (1.5) admit the entire range of L 2 -critical and supercritical 
exponents, 1 + < mi < m2 < 



Relations (1.4) and (1.5) are required to complete bootstrap estimates of the LP and 



L q norms, respectively. The arithmetic consequences of (1.4) and (1.5) are discussed 
in Section [2/71 



4. There is an adequate dispersive estimate, uniformly in to: 



sup ||e 



\Lr<t-^U\\ Lr 



(1.6) 
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where V c is the continuous spectral projection with respect to the linearized operator 



JL (depends on oj). See (2.10) below for details. We only require estimate (1.6) for 



r = p,q, which excludes endpoint cases. 



For Vq = 0, (1.6) is a consequence of the spectral assumption. For N > 2, see 
Cuccagna |11[ Corollary 2.2], built on the earlier works of Yajima [411 142] . For N = 2, 
see Cuccagna and Tarulli [T3] . 



The class of non-trivial potentials Vq for which (1.6) holds is unknown. 



All assumptions are known to be true for the monic cubic focusing equation in three 
dimensions. The lack of embedded eigenvalues is due to a numerically assisted proof by 
Marzuola and Simpson |27j . All assumptions are expected to be true for any monic L 2 - 
supercritical and energy-subcritical equation, and for their perturbations. Therefore, after 
rescaling, the assumptions should hold for sums of two monic nonlinearities in certain ranges 
of oj. Indeed, the spectral assumption is partially known for the cubic-quintic nonlinearity 
in exactly such a situation. See Asad and Simpson [2]. See Hundertmark and Lee |19j 
regarding the decay of Y+, YL . 

Theorem 1.1. Let 1 < N < oo and suppose the above assumptions are valid for uj in 
an interval I with uniform estimates. There exist positive constants oq = Oo(/, X) and 
C = C(f,X) such that the following holds. 

For any ao G (0,Oo], consider initial data no for which there exist 8q G R and ujq £ X 
with dist{ujQ,dZ) > Cao, such that: 

< a = \\u - (f)ojoe l9o \\mr\Li < «o- 

Either: 

1. (escape case) there exist u + G X, |(Jo — w+| < Cao, 0+ £ [0, 27r] and finite time 
T e xit > 0, such that, 

\\u(T e xit) - <t>u +e ld+ \\ L 2 = inf \\u(T exit ) - </> w e i9 || L 2 > 2a , or, 

2. (convergence case) there exist uj + G X, — uj + \ < Ca, t \-t 9(t) continuous and 
u = u(/) > ; such that, 

\\u(t) - ^ + e im \\ L r> < Ca(t)-t l for all t > 0. 

Recall the standard notation (t) = v 1 + 1 2 for t G K, and that p = 771,2 + 1 is the largest 
exponent in the potential energy. Since LP L 2 oc , the convergence case is an asymptotic 
stability result. Let us emphasize that the estimates must be uniform over X. In particular, 
the real eigenvalues e+,e_ and {d^^^J) are uniformly bounded away from zero. The 
restriction dist(wo> <9X) > Cao ensures that both u + G X and uj(t) G X, where u:(t) will be 



defined by (2.12). 



1.3 Context and Importance 

Asymptotic stability of orbitally stable solitary waves is well studied and has a vast, growing 
literature, initiated by Soffer-Weinstein [361 E3| an d Buslaev-Perelman [6] . 
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For unstable solitary waves, the general result of Shatah-Strauss |35j exhibits solutions 
which are initially arbitrarily close to the solitary waves but leave their neighbourhood 
in finite times. Also see Comech-Pelinovsky |10j . who give a similar result when loq is 
the borderline between stable and unstable branches. There are also results showing the 
existence of stable manifolds, solutions which converge to the unstable solitary wave, see 
[391 [261 [33] . Solutions on a stable manifold are necessarily nongeneric. Indeed, there are 
few results addressing all solutions with initial data in a neighbourhood of unstable solitary 
waves. There are three exceptions: 

1. Small solitary waves obtained from a linear potential. 

A complete description of asymptotic behaviour is known in some cases. See [38, 32J. 

2. Solitary waves of the pure-power L 2 -critical nonlinearity, g(u) = \u\nu. 

There exist disjoint open subsets K^ 1 C H 1 such that the solitary waves belong to 
JC + n Kr . In some cases, we know that solutions in /C + scatter. See Killip, Tao and 
Visan |22j and papers that refer to it. On the other hand, solutions in K~ blowup in 
finite time, as proved by Merle and Raphael [29J following a couple decades of careful 
asymptotic arguments. These blowup solutions are precisely described in terms of 
a soliton profile and a tracking error. The tracking error is arbitrarily small, and 
converges, in L 2 . Indeed, it converges in H 1 outside any ball of fixed radius around 
the blowup point. The primary growth of H 1 norm is captured by the soliton profile 
for which 11^)11^1 -> oo as u(t) -)• oo. 

The L 2 -critical nonlinearity is a degenerate case with physical relevance. The related 
literature is very large. 

• 1 4 

3. Solitary waves of the ffz-critical nonlinearity, g(u) = \u\ N ~ 1 u. 

There exist disjoint open subsets JC^ C H 1 such that the solitary waves belong to 
/C+ D K,~ . Moreover, solutions in /C + scatter as t — > oo and radially symmetric 
solutions in K," blowup in finite time. See the combined works of Duyckaerts, Holmer 
and Roudenko [181 El IS] , expanding on the energy-critical argument of Kenig and 
Merle 115]. For general energy-subcritical pure-power nonlinearities, see Akahori 
and Nawa p]. 

The nature of the blowup dynamic is unclear. Numerical results suggest a stable 
(almost) self-similar blowup regime. It is then natural to ask if there are solutions 
that blowup by following the soliton manifold, as in the L 2 -critical case. 

For slightly L -supercritical nonlinearities, g(u) = \u\ N u, Merle, Raphael and Szef- 
tel [ 3"T] have shown that the L 2 -critical blowup regime survives as sets of initial data 0, 
open in H , for which blowup occurs and can be described in terms of a member of 
the soliton family and a tracking error. At blowup time, the tracking error converges 
in all subcriticaQiorms to a fixed residue which is outside the critical space. This 
agrees with the more general result of Merle and Raphael [30J that the critical norm 
of radially symmetric blowup solutions is unbounded. Moreover, there is a universal 
lower bound for the size of the residue in L 2 oc . Should blowup with a soliton profile 
occur in any other L 2 -supercritical problems, Theorem |l . l| suggests there should be a 
similar L 2 oc lower bound on any residue. 



1 For a pure-power nonlinearity, g(u) = \u\ p 2 u, the critical norm is L 2 
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Should equation (1.1) lead to blowup, it may be structurally pertur bed by a vanishing 
multiple of u\u\ m2 ~ l to be globally wellposed. As a result, Theorem 



1.1 



shows that the 

qualitative dynamic near a branch of unstable solitons is universal, irrespective of whether 
blowup eventually occurs. We do not exclude the possibility of solutions that blowup with 
a soliton profile, following an unstable branch of solitons at some distance. Should such 
solutions exist, Theorem |1.1| suggests that the blowup dynamic is only observed once they 
lie outside a particular neighbourhood of the manifold. 

While we are concerned with large solitons, our approach will be similar to the small- 
soliton case. Kirr and Mizrak |25} [23| IM] consider the nonlinear Schrodinger equation 
with potential in dimensions two to five and detail the convergence to a center manifold of 



small stable solitons. Since they do not require (1.5), their work admits a larger range of 
nonlinearities down to the Strauss exponent. The technique in all three papers is focused 
on time-dependent linear operators, which we avoid, and is strictly limited to small solitons. 
Recent work of Beceanu [3j may offer a new route to soliton stability results in the in- 



setting. We do not know if Theorem 1.1 is optimal. 

Our approach recovers the asymptotic stability of large solitary waves with no non-zero 
eigenvalue, for radial perturbations: 

Proposition 1.2. Let 1 < N < oo, suppose that for a; € X the linearized operator JL has 
eigenvalue zero with multiplicity two, no other discrete spectra, and that all other assump- 
tions are valid with uniform estimates. Then the result of Theorem holds. Moreover, 
only the convergence case occurs. 



For small solitons, potentials Vq with linearized operators as in Proposition 1.2 exist. 



Proposition 1.2 then recovers a range of nonlinearities covered by Kirr, Mizrak and Zarnescu. 



2 Decomposition and Algebraic Relations 

First, in Section |2.1[ we properly introduce the linearized operator JL. Second, in Sec- 



tion 



2.2 we decompose the solutions u(t). This will allow us to phrase the bootstrap 



argument for Theorem |1.1[ explained in Section 2.3 We then introduce particular tools 



in preparation for the following chapters. We state the dynamic equations of the modula- 



tion parameters in Section |2,4[ the tracking-error equation in Section 2.5 and a Buslaev 



Perelman decomposition and estimate of the continuous spectral projection operator of 



JL in Section |2.6| Finally, in Section 2.7 we introduce decay-rate constants and verify 
associated arithmetic. 

2.1 Linearized Operator 

Expand the potential around (f>, 

g{<l> + e)=g{4>) + Dg{<t>)E+N(e) 

= cf>f(\<p\ 2 ) + (V(H 2 ) + 0/'(|0| 2 )(# + fc)) + M(e). (2 ' 1} 

Immediately, we recognize the linearized operator 

C(e) = -Ae + V s + loe - (e/(|0| 2 ) + 4>f'M 2 )(^ + 0e)) 
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around <fi, which satisfies idt£ = C(e) + 0(e 2 ). Equation (1.2) allows us to estimate the 
nonlinear term (e.g. \24\ equation (22)]), 



\Af(e)\ < A, |</>r~ 2 \e\ 2 + A 2 \4>\ m *~ 2 \e\ 2 + | e | mi + \e\ m * 



(2.2) 



where Aj = if mj < 2. 

Remark 2.1 (Complex- valued Functions as Vectors). Consider complex-valued functions 
u = m + iu 2 and v = v\ + iv 2 , which we write as M 2 -valued functions: [full = Ul 

[u 2 _ 

Vl . The correct inner product is 

V2 



and 



r 



Ul 




Vl 


U 2_ 




v 2 



(ui,vi) + (u 2 ,v 2 ) 



(uivi + u 2 v 2 ) dx. 



We include the complex conjugate of v since we will later consider C 2 -valued functions. We 
denote the symplectic operator by 

" +1 
-1 



Write e = E\ + ie 2 , and represent e as [e] 



. Then, 



~L+ " 




ei 


L_ 




£ 2_ 



L[e] = ((-A + W )I + V) [el 



where 



L+ = -A + Vo+ui — (/(0 2 ) + 20 2 /'(0 2 )) , L_ = -A + V + u - f{4> 2 ). 
The term I denotes the identity, and the potential term V has localized support and is 



written separately in anticipation of equation (2.19). Note that L is self-adjoint, and that 
the linearization of ( 1.1 ) near will feature JL. The kernel of JL can be found by inspection: 
L-d> = and L+d^d) = —<b. In vector notation: 



JL 



du4> 




"cf 


, JL 


"0" 







A 







0. 



Let Y + = Y re + iYim be the eigenfunction associated with e+: 

jL[y + ] = e+ [y + ] = 



-e+Y ir 



(2.3) 



Let Y- be an eigenfunction associated with e_. From (2.3), one can verify Y- = Y+ 
and e_ 



-e + . Since L_ is non-negative, and L_y, m = e+Y re , we note that (Y re ,Yj, 
Yim) > 0. Without loss of generali 
the discrete eigenfunctions of JL. Note that 



(L_Yi m , Yi m ) > 0. Without loss of generality we specify 2 (Y re , Yi m ) = +1. This describes 



P±],J 







and 







0. 



(2.4) 



S 



2.2 Orthogonality Conditions 

We assume the decomposition 



(2.5) 



where the modulation parameters u)(t), 0(t), a(t), b+(t), b-(t) 6 M. are continuous functions 
determined by enforcing orthogonality conditions^] To determine the parameters: 



a 



<W,Jp+l) 



(2.6) 

(2.7) 

(2.8) 
(2.9) 



Parameters and a(t) are not independent; see the proof of Lemma 2.2, below. We 
choose to fix a(t) = 0. As a consequence, the linearized operator JL, the eigenfunctions Y + , 
Y-, and their associated eigenvalues, are all themselves functions of time through u)(t). To 



simplify notation, this dependence is usually omitted, as in (2.5). When we consider a fixed 



operator, associated with some fixed value uj(T), we will refer to the associated linearized 
operator as and the eigenvalues as e± )W m. 

The chosen orthogonality conditions are with the eigenfunctions of the adjoint of JL. 



These conditions will allow an easy derivation of the dynamical equations in Section 2.4 
Indeed, r] is the projection onto the continuous spectrum of JL, 



r} = V c 



-i6(t) 



u(t) 



(2.10) 



whenever the orthogonality conditions uniquely determine 77. 

To see that this is the case, let Xq denote the JL-invariant subspace associated with 
the generalized kernel, and X\, X c the subspaces associated with eigenvalues e± and the 
continuous spectrum respectively, so that L 2 = Xq © X\ © X c . Let Vq, V\ and V c = 
Id — Vq — V\ denote the projection operators onto Xq, Xi and X c respectively. Explicitly, 



Pi in 



1 



2 UuJ 



2 

L- 



(UUY+l) 



I/], J 















(UUY-j) 



(2.11) 



HY+i, 



2 {Y re , Yiu^j 2 (Y r( 

and we note that V c \ f\ satisfies the orthogonality conditions. 

Lemma 2.2 (Ability to Modulate). Fix uj\ £ I, 9\ £ M and Banach space L r , 1 < r < 00. 
Then, foru in some neighbourhood of (f)^ 1 e idl , there exists a Lipschitz map a4 (9, oj, &_)_, 6_) 



such that, under decomposition (2.5), orthogonality conditions (2.6)-(2.9) are satisfied. The 
radius of the neighbourhood is of the order dist(wi, dl) as dist(wi, dT) — > + . 



2 By way of comparison with Buslaev-Sulem, note that the coefficients b±(t) and the components of [Y±] 
are real- valued. In the case of stability and eigenvalues ±iA, the second component of the eigenfunction in 
vector form is purely imaginary (corresponding to an entirely real- valued eigenfunction) , and the appropriate 
decomposition is z(t) ftp} + z(t) [^] . 
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For our application we take r = p = + 1. In particular, |u>o — w (0)| < a and the 
modulation parameters are continuous in time. 

Proof. Define the map p:(9,w,u)elxIxL r 4 p(6, u, u) G R 2 by 



p(6,u>,u) 



lei J 







WTO 





where e = ue 



-w 



At (6*i , wi , 0oji e 101 ) , /0 = 0, and the Jacobian with respect to (9,ui) is 



2^a 



L-l 



|^w||VllL2 |(( . : , 




The Jacobian is nonsingular and its inverse is uniformly bounded over uj\ € X by assumption. 
We may apply the implicit function theorem on Banach spaces (e.g. Berger |5j) to solve 



for (8(u),uj(u)) such that p(6(u),u(u),u) = 0. We then define rj, b+, and 6_ by (2.10) and 

<27TTb. □ 



2.3 Critical Time T crit & Proof Strategy 

For all data under consideration, there exists Td ec > such that we may decompose the 



solution in the manner of (2.5) for t G [0,T<2 ec ), 



u{t) = (0 w(t) + b+(t)Y+ + &_(t)y_ + r,(t)) e ie ®. (2.12) 
Define a new time scale T cr it < r^ ec in terms of persistent good control of 6+(t), 

T crji = sup {ViG[0,T], |6+(t)| < a(t)" 1 }. (2.13) 

T<T dec 

Our proof proceeds on two paths, depending on T cr i t : 

1. T C rit — -^dec 

This case is considered in Section [3j We prove that T cr it = +oo and that the solution 
converges to the soliton family. 

2. Tcrit ^ Tdec 

This case is considered in Section |4j We prove that the growth of b+(t) cannot be 
controlled, and that exit from the neighbourhood of the soliton family occurs at T ex n, 

with T cr it < T ex it < Tdec- 

2.4 Dynamic Equations 

Substitute ( |2.12 | into (1.1) and use the expansion (2.1): 

i(6 - u) ((f) + e) + Cod^ + d t e = -i (Ce + M{e)) , 



where e = b+Y+ + b-Y- + rj, and M was defined in (2.1 ). In vector form, 



o 



+ (e-u) 



(2.14) 



JL Ifel + 



u)J[e}+JN{e), 
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where we use the notation N(e) = [A/"(e)J and write dt [e] in full to emphasize its terms 

Ol 



will be handled differently. Take the product of (2.14) by J 




and use (2.6) to integrate 



dtMJ 



by parts in time. We get 



to (d w <j),(j)) =Co (M 
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+ N{e) 



(2.15) 



Similarly, with —J 







and using (|2.7|), we get 




2 







N(e) 







(2.16) 



Now with J [1+] and using (2.8) to remove the terms in dtr\ and Lt/, we obtain (recall 
that 2{Y re ,Y im ) = 1): 



6_ =e-b. 



+ 



'u}({fnj,d w 3 [Y+j) - (b+d^ {Y+j + b-d u lY-l J (2.17) 



2(Y re ,Y im ) \ +([0-u})[e}+N(e),\Y + } 



Finally, with J [YL] and using (2.9) to remove the terms in dtr\ and L77, we get: 
b+ =e + b + 

1 /> ((fo]|,0 w J [Y_J) - (6+0* [Y + ] + 6_0 W [YLJ, J [Y_]))^ (2 . | N) 
" 2(Y re ,Y im ) I + / (0 - ^ [ e ] + N(e), [Y_] 



2.5 Tracking-Error Equation 

It will prove more straightforward to estimate the tracking-error 77 in terms of a fixed 
operator. Let V c t denote the projection onto the continuous spectrum of the operator with 
some fixed w(T)j^] Denote rj = V Cs Tf}- We do not change our choice of decomposition. 
We first isolate linear terms in 77 in (2.14), so that 



d t [rjl = JL [77] + 0- to) J [77J + B + JN, where 







+ i JL - d t + 



^)j}(6+ [Y+l + 6- [y_]). 



Note that V = L — (—A + Vq + w)I is a localized potential. We further rewrite 

d t [77] = JL T M + (0 - WT )J [r?] + J(V - V T ) M + £0 + JiV. 
Applying "P^y to all terms, we get 

a [tTI = JL T {rj\ + (O - lot) Vc,tJ H 

+ V c ,t (- (e - ut) J (V c ,t - V c ) [77] + J(V - V r ) [77] + B + JN 



(2.19) 



The convergence of cu(t) as t — > T cr u will be established by Lemma 3.1 
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where V c ,tB = V c ,tBq and 
B = - (TV - Vc 







+ 



UJ 







(V C , T - Vc) {(b + - e + b + ) {Y+j + (6_ - e_6_) [K.]^ 

w (b+d u {Y+j + [y_]) + (e - J) J (6+ + &_ [r_ 



Observe from (2.11 ) that both V — Vy and Pc.t — V c are localized potentials, of the order 
\ui — ojt\, depending on I and provided \co — ujt\ is sufficiently small. Also note that B has 
localized spatial support, and 



B = - ur\ + |6+| + |6_|)(|w| + 
2.6 Buslaev-Perelman Estimate 



CO 



+ 



e+6 H 



+ 



e_6_ 



(2.20) 



Later we will assume the tracking-error rj is small, and then the evolution given by (2.19) 
should be essentially linear. To handle the leading order correction in (0 — u), we use a 
technique introduced by Buslaev and Perelman |6j , later discussed by Buslaev and Sulem [7] . 
For N >3, the proof is given by Cuccagna |12| . For N = 2, the proof is claimed by Cuccagna 
and Tarulli [13]. Let V+ and V- denote the spectral projection operators onto the positive 
and negative continuous spectrum. That is, V C ,T = V+ + V- ■ 

Proposition 2.3. 

- p CjT J = P CiT J _1 = iV+ - iV- + K, (2.21) 
where K is a localizing operator, bounded L q \-t L r , for any pair l<r'<2<q<oo. 



Remark 2.4. To motivate (2.21), note that an equation of the form c\ [/] = JL-p [/] may 
also be written in the form 

idt 

where the potentials Vj are due to both Vo and the nonlinearity. We may view this matrix 

A + uj + Vq 
A-u-V 



7" 




'-A + u-Vi V 2 




7" 


/ 




V 3 A - uj + V A 




/ 



, for which (2.21) is an identity 



operator as a perturbation of 
when K = 0. 

Let us decompose [77] = V C: t according to V± , and incorporate the accumulated 
error in tracking the phase, 



V± = (V± 



The evolution of rj± follows from (2.19) 
d t v± = JL T ?j ± + V±e ±ie ^ 



- uj T K [rjl + J(V - V T ) M + B + JN ( e ) 



(2.22) 



where we have abused notation to absorb (6— lot) CPc,t — Vc) i n t° V — V^. Equation ( 2.22[ ) 
will be our means to establish L p and L q estimates for rj through the equivalent Duhamel 
formulation: 



V± 



JJLn 



3 (t-s)JL TV±e ± t e(s) J^_ w j K [r|+J(V-V r )i+B + JJV (e)) (s) ds. (2.23) 
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2.7 Arithmetic 



For later reference, we collect here an assortment of arithmetic facts. Due to ( 1.4 ) and ( 1.5 ), 
there exists 5 > sufficiently small such that we may define: 



for Or, > 



m + 1 

2 2 

1- 5 for Up < , 

. mo + 1 mo + 1 

AT 

with cr g < min{l, ~ 1) ~~ !}■ 



(2.24) 



Note that Ug G [<r p , 1), corresponding to some g E [p,p m ax)- Such a choice is only possible 
for N > 2, since 



(mo + l)o- g > 2, 



(2.25) 



and, in particular, a q > \. Constant 5 > is included to ensure strict inequality for (2.25), 
and is otherwise neglected by taking it sufficiently small. Let us emphasize that 



mi > 1 + — (1 + a q ). 



(2.26) 



Consider m^B^o r , 
Then 



N 



mj — 1) — <7 p , an expression that will arise during norm interpolation. 



mj8j(jp > 1 and TrijOj > 1, 



(2.27) 



which are due to (1.4) and a p < 1, respectively. Due to q < p ma , x , 



2 ( I - 1 ) < l + ~ < 



(2.28) 



Consider mjdjG p = y (mj — 1) — a q , another expression that will arise during norm inter- 
polation. Due to (2.26), and since <r p < 1 < mo, 



mj9jcr p > 1 and mjQj > 



1 

m ' 



(2.29) 



Finally, one can verify that i (l — ^) > 1 — m ■ ' "-- 



l J V 2 



j is a consequence of >m\ > 1, 



and since mf — 3mi + 2 > when mo = m\ < 2. The second expression of (2.29) implies: 



(l-%)^(^)+m,-%>l. 



(2.30) 



3 Convergence Case 

In this section, we will consider the case of T cr n = T^ec- We will prove that T cr u = +oo and 
convergence to a soliton by means of a bootstrap argument. The estimates shown below 
will be reused in Section |4j 
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3.1 Hypotheses and First Estimates 



For t G [0, Tdec), we decompose u(t) as in (2.12). Assume that Th yp G (0, T cr it] is the last 
time for which the following hypotheses hold for all t G [0, T^jm): 



|6_(t)| < 2Cia(t) -1 , 
\ LP < 2C 2 a(t)-^, 
\lo < 2C 3 a(t)- a «. 



(3.1) 



Universal constants Cj will be determined later. Recall that from (2.13), the definition of 
Tcrit, we also have 

\b+(t)\ < a(t)-\ (3.2) 



Under these hypotheses, the goal is now to obtain the same inequalities as (3.1 ) without 



the factor 2, by choosing otQ small enough. We now fix T = Th yp . First, we prove estimates 
on the parameters 6 and u. 



Lemma 3.1. For all t G [0, T% p ), we have u(t) G X, 

\uj(t)\ + \0(t) - u(t)\ < Ca m °{t)- moCT \ 

\ui(t)-w(T)\ < Ca mo (t}- moa " +1 , 
\e ±M t)-e ±MT) \<Ca m »(tr m ^ +1 . 

By exactly the same proof, we have 



which implies 
and 



(3.3) 



b± - e±b± 



< Ca mo {t) 



-m a q 



Proof. We will only prove the first inequality. To do this, we add (2.15) and (2.16), to 



obtain |w| + |# — uj\ < (\6j\ + \8 — uj\)0{e) + |(7V(e), potential) |. With ao small enough, we get 
\ui\ + \6 — uj\ < (N(e), potential), and (3.3) follows from ( |2.2[ ) and the bootstrap inequalities 
(3.1) and (3.2). We use the L q control of rj since it may be stronger by (2.24). □ 



3.2 Improved Estimates 

Estimate on 6_. We can now improve the estimate on Fix e\ = e_ jU1 (T) < 0- From 
Lemma 13. II we obtain: 



e\b- 



< a mo (t)- mo ^ ((t) \b-(t)\ + 1) < Coa mo (t)- moa i, 



where Co > is an explicit constant depending on X. Moreover, since ||£(0)||#i n £,i = a, we 
may assume |6-(0)| < C§a. Integrating between and t, 

|6_(t)| <e +eit C a+ f e +eiit - s) C a mo (s)- moa " ds 
Jo 



< ae +eit 



+ a 



mi) 



1=1 



e +e.(t-s) ds+ f* ( S )- m ^ds ) , 



with the obvious corrections to integral bounds when t ~ 0. For oq sufficient small, we 
have shown |6-(t)| < a (1 + ln(i)) (t)~ m ° aq < C\a{t)~ l where C\ > Co is the appropriate 
universal constant. 
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Decomposition ofrj. Recall from (2.10) that rj = V c r), implying 



■q = rj + {Vc- V c ,t) V, 



so that norms of rj, rj, and rj± are all comparable. Due to (2.6), we have 



\u(t)\\ 2 L 2 « (4>oj + -n,4>uj + 77) 



|2 1 11 ||2 



where we have ignored terms of order b±(t). From Lemma 3.1 we conclude 

m 

| L 2 < C 4 a^, 



(3.4) 



where the constant C4 depends on X. Now we estimate terms of (2.23), in some norm U , 



with r' G (1)2). From Proposition 2.3, (2.25), and Lemma 3.1 

-oj t ) KM \\ Lt/ < (e-ur) U\l* < Ca 2+ (t)^ mo+1 ^ +1 < a 2+ (t)^ 1+ \ (3.5) 
Note that we use here standard notation 1+ to designate a number slightly bigger than 1. 



The same estimate holds for both ||J(V — Vy) [77] \\ Lr ' and ||-B|| ir /. Indeed, from (2.20), we 
have 

\\B\\ Lr > < Ca 2+ (^t)- 2mo ^ +1 + ( t )-»»o»«-i) ; ^3 g) 

which is lower order provided a q < 1. 



LP estimate ofrj. To start, we estimate ||-/V(e)|| iP / using equation (2.2) 



\lp' 



<(\b + \ + \b\ 



+ \\v\ 



which may be interpolated as 



1 



1 



\ri\\ L L jP > < WvWh WvWJ : 



1 1 



(3.7) 



with ; = — H or, equivalently, 9j = — n h^^ We note that 0,- G (0, 11 provided 

m jP ' p 2 1 - I 

1" 



2 J 1 J < rrij < p — 1, which is implied by 1 + < mj < 777,2. With (3.1), (3.2) and 



(2.27), we have 



If rrij > 2, 



^\\ LP ><a 1+ (tr {1+) - 



\A 3 r^ 2 e 2 \\ Lv > < (\b+\ + |fe_|) 2 + \\ V \\h < a 2 (t)- {1+) , 



(3.8) 



which we view as a correction to (3 



Apply (1.6) to all terms of (2.23), with r = p. Since L p is energy subcritical, we may 



use that ||t7±(0)||#i is small to improve the bound on the linear term, 



„<JL r ~ 



Co 



V±(0)\\l P < ~fa{t)-^. 



The universal constant C2 is determined by this relation. For the other terms, apply (3.5), 

< ^a{t)-^ + Ca l+ J* \t - s\- ff >(s 



( 3.6[ ) and to get 

\\v±V)\\lp 



-(1+) 



ds. 



(3.9) 
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This proves ||??(i)||ip < C2a(t) ap , as desired, by assuming oq > is sufficiently small. 



L q estimate ofrj. We consider now the L q norm of r\. From (2.23), we have 

¥\±{t)\W < l|e UL ^ ± (0)|U, 

+ J \\e^ JLT V ± ((# - w) K [77I + J(V - V T ) M + S) \\ Lq ds 



+ 



[ l \\e {t - s)JLT V±N(e)\\ Lq ds 
Jo 



I + 11 + III. 



(3.10) 



To estimate these terms, we use dispersive estimate (1.6). We first get 



(3.11) 



The universal constant C3 is determined by this relation. 

Term III is treated in a similar way as for the LP estimate. We have 



|^||^<(|6 + | + |6-ir + ||ry| 



where now we interpolate according to 



3 < 

\ L va ]q l - 



\'l\\i2 W'IWlp i 



(3.12) 



with ' . = 1 " J ' ""' 



1 1 



rrijq 
1 



+ — , or, equivalently, 9j = — \y e no t e that E (0, 1] provided 
2 P 2 ~ p 

1 



2 1 < rrij <p[ 1 . The former inequality is (2.28), while the latter inequality 

is true for any q > p > nij + 1. From (3.4), the terms of (3.12) give a with the exponent 



(1 — 6j)rrij^Y' + m i^ji which is greater than 1 by (2.30). We have 



\e m i\\ Lg ><a 1+ (t)- m ^, 



(3.13) 



and if mj > 2, 



lA^-VH^ < (|6+| + |6_|) 2 + hill, < a 2 (t)- 2 ^, 



(3.14) 



and both decay faster than (t) , due to (2.29) and (2.25), respectively. With (1.6) 



III < a 1+ f \t - s\-°* {s)- {1+) ds < a 1+ (t)- a ". 
Jo 



(3.15) 



The integral is a simple calculation since a q < 1. Term II may be included with III. Its 



contribution is controlled by (3.5) and (3.6) applied with r = q 



Assuming ao > is sufficiently small, we have shown that 



-a <C 3 a(t)-^. 
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3.3 Bootstrap Conclusions 

It was shown in the previous section that, for all t G [0,Th yp ), 

IM*)I < c ia (t)-\ 

\\v(t)h* < C 2 a(t)-°>, 
\\ri(t)h« < C 3 a{t)-°<. 



(3.16) 



As these are continuously evolving quantities, by (3.1 ) it must be that Th yp = T cr i t . Together 
with the estimate for 6 + , this proves proximity to the soliton family as t — > T cr n, and so 

T crit . Therefore, if T crit 



by Lemma 2.2 the decomposition can be extended past t 
then Td e c = +oo. 



Tde 



By integration of the first equation of (3.3), there exists oj + G X, |w(0) — w+| < a m ° such 



that the remaining estimates of (3.3) hold. It is essential here that m§o q > 1. From (3.16), 



we have that \\u(t) — 4>uj(t) e \\lp < Ca(t) fTp . Together we have, for all t > 0, 



je(t) | 



\\u(t) - <j>u + 

For Theorem [lTJ we use fi = min{<T p , m§o q — 1} > 0. 

4 Escape Case 

In this section, we consider the case of T cr it < T^ ec . The arguments of Section [3] apply for t G 
[0, T cr it). We extend these arguments to prove that parameter |& +>aj (()| grows exponentially 
for an interval of time after T cr u- 

4.1 Hypotheses and First Estimates 



Assume that T ex u G (7c r ^,Td ec ] is the last time for which the following hypotheses (4.1) 



(43| hold for all t G [0,T exit ): 

\b-(t)\ 

Ht)\\LP 

\\v(t)h, 



< 2C 5 |6+(t)|, 

< 2C 6 {\b + (t)\ mo + a(t)-^) 

< 2C 7 (|&+(t)| +a(t)- ff «). 



(4.1) 



Universal constants Cj will be determined later. In place of (2.13), we make two further 
hypotheses: for all s < t G [T crit , T exit ), 

b+(t) 



e 2(t-s) ^ 



b + (s) 



< e2 



a e 2 (t~s) 



where e 2 

e -MTcrit) 



e +,w(T crit )i the positive eigenvalue of for T = T cr it- 
—&2 < 0. Our final hypothesis is, for all t G [0,T ex u), 

\b+(t)\ < 3a . 



(4.2) 

Recall that e± = 
(4.3) 



As a particular consequence of (4.2) and (2.13), note that for t G (T cr it,T ex it), 

\b+{t)\ > \b + (T mt )\ > aiTcrit}' 1 > a(t)-\ (4.4) 
Under these hypotheses, our goal is to obtain hypotheses (14. 1|) without the factor 2, and 



hypothesis (4.2) with tighter exponents. We will argue that hypothesis (4.3) cannot be 



improved, and conclude the expected exit behavior. 
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Lemma 4.1. For all s < t E [T cr it, T ex n), ui{t) £ X and: 



|w(f)| + |0(t) - < C (\b+(t)\ ma + a mo (i)- moff9 ) , and 

\co(t) - lo(s)\ < C (|M*)P + a mo (s)- mo<J?+1 ) , which implies 
\e± M t) ~ e ±Ms) \ < C (\b + (t)r + a m Hs)- m ^ 



Proof. The argument is the same as for Lemma 3.1, based on the following estimate: 
| (N(e), potential) | < |6+(t)| mo + a m ° {t)~ moe7q . When integrating quantities in \b+(t)\ r , for 
any exponent r > 0, we use (14.21): 



(4.5) 
□ 



|6 + (r)| r dr < \b+(t)\ r / e-i re2( *- r) dr. 



4.2 Improved Estimates 



Growth estimate for 6+ . By Lemma 4.1 the dominant forcing term of (2.18) is (N(e), [fK-J), 



and it reads: b + — e2b + = O (|6 + | m ° + a m °(t) mocr <^_ Taking ao to be sufficiently small, we 
may assume — e2^+| < |e2|6+|. After integration, this is a stronger statement than (4.2 ): 



b+(t) a. 
b + (s) 



e fe 2 (t- S ) < »+W < e |e 2 (t- S )_ 



(4.6) 



Estimate on 6_. The same argument applied to (2.17) gives 

\b-(t)\ < e+^- T ^) |6_(T crit )| + f e + e i(*-) |6 + (r)| dr. 



Compare the estimate for 6_ from the previous chapter (3.1) with (2.13) to see that 



\b-(T cr it)\ < C\ \b + (T cr it)\. We conclude that |6_(t)| < C5 |&+(t)| for some C5. 



Decomposition of r\. Analogous to (3.4), and using ( 2.4 )-rt2.9|), we have 



\\v(t)\\ L 2<C 8 {\b + (t)\^ +a^). 



(4.7) 



Under our new hypotheses, let us revisit (3.5 ) and (3.6). As before, we use (2.20 ), Lemma 4.1 



and (2.25) to conclude: 



LOT ) K 



+ ||J(V-V T ) 



+ \\B\\ Lr , < [\b+{t)\ mo + a mo (t)- (1+) 



■ (4-8) 



LP estimate of n. For t G [T cr it,T ex i t ), we estimate ||iV(e)|| iP ' using (3.7) and (4.1): 



\\N(e)\\ L r> z \b+r + £ iMi£"' i)B *cr'' 0M*)r + a<r CTp n . (4-9) 



Due to ( |2.27D and (|L3]), we may bound ([49]) by ( |fc+| m ° + a 1+ (ip {1+) ), for any universal 



constant C9 > 0, by taking the universal constant Cq sufficiently large and ao sufficiently 
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small. Now apply (1.6) and use both (4.9) and, for t G [0,T cr i t ), the estimates from the 



previous chapter that led to (3.9): 



\v±(t)\\ L p < ^a{t)-^+Ca l+ f \t - s\~ a " (sr (1+) ds 
* Jo 

+ 7^ I \t-s\-°*\b+{s)\ mo ds. 

°9 JT crlt 



From (3.9), the first terms are bounded by C2a{t) ap . Assume that > Ci. Integrate the 



final term with (4.5 ), and take Cg to be the constant factor. This completely determines Cq. 



We have shown \\fj±(t)\\ LP < C 6 (a(i) _<Tp + \b+(t)\ mo ). 



L q estimate ofrj. Recall the terms I, II and III of (|3.10). The estimate (3.11 ) of term I 



still applies. For term III, we first consider (3.12) for t > T cr u, using (4.1), (4.4) and (4.7), 
IMI^y < o(l) (\b + r + a(t)-^r dl < \b + (t)\ 1+ , 



where the second inequality relied on (2.29) and ( |2.30 ) for the exponents of b + and a, 
respectively. If rrij > 2, 

WM^-^Wi* <\b + \ 2 + a 2 (tr 2 ^, 



as in (3.14). As before, all terms decay faster than (t) and we have: 



III< 



f \b + (s)\ 1+ ds + a 1+ f \t-s\-^{s) 



-(i+) 



ds, 



which are the same integrals as (4.5) and (3.15), respectively. Term II may be included for 



s G [0, t), controlled by a combination of (3.5), (3.6) and (4.8): 

11 + III < \b + (t)\ 1+ + a 1+ (t)-^. 
This concludes our estimate of ||^±||i,9- 

4.3 Bootstrap Conclusions 



In the previous section, we proved (4.6) and that, for all t £ [T cr n, T ex n) 



\b-(t)\ < C 5 \b+(t)\, 
\W)h* < C 7 (\b + (t)\ + a(t>-'«) , 
\ LP < C 6 {\b + (t)r+a(t)-°*). 



(4.10) 



These are continuously evolving quantities. From (4.3), (4.4) and d4.6ft, we conclude that 



Texit < T crit + - — In (3—(T crit ) ) < +oo. 
4e2 V a 



Assuming oq is sufficiently small, Lemma 2.2 (4.3) and (4.10) imply that the decomposition 



can be extended past time T ex n, and hence T ex n < T^ ec . The only possible failure at time 
t = Texit is (4.3), and so we conclude that b+(T ex u) = 3ao- The conclusion of Theorem |1.1| 
then follows for some w+ ~ co(T ex n). 
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